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Abstract: We present a method to correct for deflections of ultra-high energy 
cosmic rays in the galactic magnetic field. We perform these corrections by 
simulating the expected arrival directions of protons using a parameterization 
of the field derived from Faraday rotation and synchrotron emission measure¬ 
ments. To evaluate the method we introduce a simulated astrophysical scenario 
and two observables designed for testing cosmic ray deflections. We show that 
protons can be identified by taking advantage of the galactic magnetic field pat¬ 
tern. Consequently, cosmic ray deflection in the galactic field can be verified 
experimentally. The method also enables searches for directional correlations 
of cosmic rays with source candidates. 



1 Introduction 


The magnetic field of our galaxy is presumed to cause substantial deflections of ultra- 
high energy cosmic rays prior to their observation. To understand the deflections, 
the directional dependence of the field pattern and the directional variations in the 
magnitude of the field are most relevant. 

Explicit predictions for deflections depending on the cosmic ray arrival direction, 
its energy and charge can be obtained through recent parameterizations of the shape 
and magnitude of the galactic magnetic field [D El El la E]- They are based on 40,000 
Faraday rotation and synchrotron emission measurements. Exploring cosmic ray 
deflections thus constitutes an independent test of the field parameterizations, and 
can in addition be related to obtaining information on the origin of cosmic rays. 

In this contribution we present an analysis method to investigate cosmic ray 
deflections in magnetic fields. The method is based on propagating protons through 
the field to obtain their expected arrival directions. Owing to the forward simulation 
technique we can test the validity of a field parameterization without imposing 
potentially unphysical conditions on measured cosmic rays. 

Our analysis concept consists of testing a combined model with three assump¬ 
tions: The first assumption is that coincident directions of cosmic rays and source 
candidates exist outside of our galaxy, before the cosmic rays are deflected in the 
galactic magnetic field. The second aspect of the model is that the above parameter¬ 
izations of the galactic magnetic field reflect the field characteristics and magnitude 
as realized in our galaxy. The third assumption is that the magnetic field is suffi¬ 
ciently strong to separate protons from nuclei through small angular deflections. 

To assess the method we simulate an astrophysical scenario that we call ‘data’ 
throughout this contribution. To generate these data we use the cosmic ray prop¬ 
agation program CRPropa version 3 [B]. As sources of cosmic rays have not been 
identified yet, we make the following choice for the simulated data: We assume ac¬ 
tive galactic nuclei (AGN) within the GZK horizon [7] to be source candidates, and 
use neutrinos as navigators for selecting source directions. 

This contribution is structured as follows. We first explain the simulated astro- 
physical scenario and the corresponding simulated data set. We then present the 
method of calculating the expected arrival directions of cosmic rays after traversing 
the magnetic field. We analyze the angular distances between the observed cosmic 
ray arrival directions and the expected arrival directions. We also quantify the prob¬ 
ability of finding several cosmic rays associated to single AGNs. We then study the 
effect of alternative model assumptions when analyzing the simulated data. Finally 
we present our conclusions. 
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2 Astrophysical scenario 

Here we assume that a subset of active galactic nuclei (AGN) can be considered as 
sources of cosmic rays. As navigators to select these AGN sources experimentally 
we use the arrival directions of 24 high energy neutrinos published by the IceGube 
Gollaboration jH]. These neutrinos have energies above E = 4Q TeV, and were not 
assigned strong evidence for resulting from atmospheric background. 

For the AGNs we use the VGV catalogue [9] and consider AGNs with distances 
2 ; < 0.018. This value is motivated by the analysis performed in [10] and provides 
sources within the GZK horizon. With the criterion of selecting the nearest AGNs 
of the 24 neutrinos we obtain 22 AGNs where two neutrinos have the same AGN. 
The AGN sources are shown by the star symbols in Fig. ig>. with the color code 
indicating the AGN distances. 

Our aim is to obtain a simulated cosmic ray data set corresponding to A^ = 231 
ultra-high energy cosmic rays with distributions similar to the one published by 
the Pierre Auger Gollaboration m- For this we combine 10% of an astrophysical 
scenario simulated with the GRPropa program jH] with 90% isotropic cosmic rays 
following the geometrical acceptance of the observatory na. 

For the 10% contribution of the astrophysical scenario we generate 10^ cosmic 
rays at each of the 22 selected AGN sources. As the initial composition we use a flat 
distribution of nuclei with charges between Z = 1,...,26. Their minimum energies 
are = 50 EeV, and their maximum energies correspond to their rigidities Z x E^ax 
with Emax = 500 EeV. 

Upon arrival at a 0.5 Mpc observer sphere around our galaxy we use galactic 
magnetic field lenses to project the arriving cosmic rays onto Earth. The lenses 
consist of matrices based on the HEALPix format [13] where we divide the sphere 
into 49,152 equally sized pixels of approximately 1 deg. The matrices contain the 
probability of a cosmic ray entering the galaxy with energy E at direction i to be 
observed in direction j. The technical details of the lenses and their production are 
outlined in na. The lenses used in this contribution have been calculated with the 
GRPropa program and the JF12 parameterization of the galactic magnetic held [2]. 

We then select cosmic rays according to the geometrical acceptance of the ob¬ 
servatory [12], and an energy distribution corresponding to the measured energy 
spectrum [15]. Of these cosmic rays we randomly select the 10% contribution of 
cosmic rays arriving from the AGN sources where we require 70% to be protons. 
The total sample of the 231 simulated cosmic rays therefore contains 7% proton 
signal and 93% background (grey circular symbols in Fig. [^). 
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3 Analysis method 

To investigate the galactic magnetic field using the cosmic ray data set described in 
the previous section, we simulate their expected arrival directions assuming that all 
cosmic rays are protons. We start these proton simulations with the same set of the 
22 AGNs mentioned above. 

To take into account effects of extragalactic magnetic fields we apply a Fisher 
probability distribution centered at each AGN direction with a concentration pa¬ 
rameter K depending on the AGN distance and the cosmic ray energy. For an AGN 
at distance D = 10 Mpc and a proton with E = 52 EeV the angular spread amounts 
to cr = I/a/k = 6 deg. 

This probability distribution is then projected onto the Earth using the corre¬ 
sponding magnetic lens described above. To obtain a single direction for which the 
arrival probability is greatest, we calculate the radius r^o containing 50% of the 
arrival probabilities. We select the pixel with the smallest radius rso and use the 
center of this pixel as the expected arrival direction of the proton. 

The procedure for calculating the expected arrival direction of protons with 
energies E = 5, 10, 20, 50 EeV is visualized exemplarily in Fig. [^. The star symbol 
denotes the initial direction outside the galaxy. Note that the direction of the lowest 
energy proton also corresponds to the deflection of an ionized Neon nucleus {Z = 10) 
with energy E = 50 EeV. This implies that - in directions with a sufficiently strong 
magnetic field - protons in the cosmic ray data can be identified to some extent by 
small angular distances to the expected arrival directions. 


4 Angular distance and clustering strength 

In Figj^ we visualize two angular distances used in this analysis. The angle a 
denotes the angular distance between the measured arrival directions of a cosmic 
ray and its nearest AGN. The angle ocmf denotes the angular distance between 
the measured cosmic ray and the expected arrival direction of a proton with an 
AGN-coincident direction outside the galaxy. 

In Figj^ we show the cumulated number of cosmic rays arriving within a distance 
of angle cTo to their nearest AGNs. The triangular symbols present the angular 
distances between the measured cosmic rays and the expected arrival directions 
for protons including the magnetic field corrections. The histogram denotes the 
uncorrected angular distances between cosmic rays and AGNs. Below a few degrees, 
more cosmic rays appear near AGN directions when including the field corrections. 
These small angular distances are consistent with our simulations of the protons. 
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which supports the above-mentioned third assumption of our model that protons in 
the cosmic ray sample can be identified to some extent by exploiting the magnetic 
field deflections. 

In Figj^ we present the change in the angular distances a — ctcMF without and 
with magnetic field corrections. At positive values the cosmic rays come closer to 
the expected arrival directions, while at negative values the cosmic rays are nearer 
to the AGNs without field corrections. In the region with small angles we find more 
cosmic rays with improved angular distances when applying the field corrections 
compared to events for which the uncorrected distance is smaller. 

We quantify the change in the angular distances by the asymmetry 

^ _ 2 A^(q > (Ugmf) — A^(q < (Ugmf) 

N{a > acMF) + ^{oL < ocmf) 

which can take on any value from —2 to 2. Using e.g. a maximum angle of Oo = 5 
deg the angular asymmetry in the data is found to be A = 0.96. This positive value 
is a measure of the overall improvement in the angular distances between cosmic 
rays and AGNs when applying the field corrections. 

We also investigate clustering of cosmic rays with AGN directions. In Figj^ we 
present the frequencies of the cluster sizes m where we count associations of cosmic 
rays with AGNs within 5 deg angular distance. This yields configurations containing 
singlet, doublet, triplet, and sextet clusters (m = 1,2, 3, 6). The triangular symbols 
represent the data when including magnetic field corrections, and the histogram 
without them. 

To quantify the observed clustering strength we use the multinomial probability 
distribution: 

P(nu-,n.n,N - N,„) = pT - P? (1 - Pi..)"-"" (2) 

ni!...n22! (A - Nhit)'. 

The value P describes the expected level of trivial clustering between the N = 231 
cosmic rays, and the 22 AGNs where the latter are distinguished by identifiers. 
AGN i has an average hit probability of pi, and was correlated with rii cosmic 
rays. The total number of cosmic rays associated with one of the AGNs is Nhu = 
Y2- rii- The remaining N — Nhu cosmic rays without AGN correlations had a no-hit 
probability of (1—Piso)- Summing the hit probabilities for the AGNs at their nominal 
arrival directions for angular distances below 5 deg gives pj^o = 5.3%. This includes 
the geometric acceptance of the Pierre Auger Observatory [T2]. When including 
the magnetic field corrections a slight energy dependence is observed where the 
probability increases on average to piso = 5.7%. 
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When applying the magnetic field corrections to the data, the logarithm of the 
multinomial probability (|^ amounts to log]^o(-PGMF) = —15.3. Without the field 
corrections the level of clustering is smaller and results in logio(-P) = —11.5. The 
change in the clustering strength thus amounts to logio(TGMF) — logio(-P) = —3.8. 

5 Analyses using alternative assumptions 

In the previous section we performed the analysis on the simulated data set using the 
correct magnetic field and the correct set of AGN sources. With 7% signal protons in 
the data and 93% background contributions, the two observables indicated overall 
improvements in terms of the angular asymmetry A = 0.96, and the clustering 
strength logio(-PGMF) — logio(-P) = —3.8. The results are visualized in Fig|^,c,d by 
the red cross. 

Below we investigate whether the same improvements can be obtained by chance. 
In a first test we apply typical values for the expected deflections in the magnetic 
field; however, we assume neither that the cosmic rays and AGNs are truly corre¬ 
lated, nor that the above pattern of the magnetic field is correct. 

We use a simulation of 10, 000 event samples with the AGN directions, isotropic 
cosmic ray arrival directions, and random patterns for the galactic magnetic field 
corrections. The geometrical acceptance of the Pierre Auger Observatory is included 
as presented in [12]. The distribution of the change in the clustering strength versus 
the angular asymmetry obtained from these simulations is shown in Fig. by the 
box symbols, and appears to be centered at zero with a slight anti-correlation in the 
two observables. In 0.05% of the events we find equal values or improvements in the 
two observables with respect to the values obtained in the data analysis. 

In a further test we use lenses produced for correcting deflections of antiprotons 
in the galactic magnetic field. These lenses reverse the galactic magnetic field for 
protons. Applying the reversed field in our data analysis instead of the correct field 
orientation, the angular asymmetry and the change in the clustering strength both 
disappear (red cross in Fig. [^; the box symbols refer again to a simulation with 
random field directions, isotropic cosmic rays, and nominal AGN directions). Thus, 
the above-mentioned second assumption on the validity of the field parameteriza- 
tions can be evaluated in comparison with analyses using random directions and 
field reversal. 

We also investigated the above-mentioned first assumption on angular corre¬ 
lations between AGNs and cosmic rays. For this test we took the nominal field 
parameterization and the nominal cosmic ray arrival directions, but 22 chance di¬ 
rections instead of the AGNs (box symbols in Fig. |^, compared to the cosmic ray 
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data analysis denoted by the red cross). Also here the chance distribntion is cen¬ 
tered at zero, and is disjnnct from the valnes obtained with the data. A similar 
distribntion is obtained when nsing isotropic cosmic rays, while keeping the nominal 
held parameterization (JF12) and the nominal AGN directions (not shown here). 
Thns, hnding improvements at the level of onr data analysis with arbitrary direc¬ 
tions is nnlikely, implying that angnlar correlations of cosmic rays and sonrces can 
be evalnated. 

We also stndied the inhnence of a potentially limited knowledge abont the sonrce 
directions. In this test we hrst performed the analysis using the correct magnetic 
held and the nominal cosmic ray data set (red cross in Fig. [^). We then varied 
the AGN positions within an angular uncertainty of 15 deg. Such variations could 
appear, e.g., when studying direct correlations of cosmic rays and cosmic neutrinos. 
The box symbols show the results of the 10, 000 variations. Even if the source 
directions are not perfectly known, an improvement in the two observables owing to 
dehections in the galactic magnetic held can be observed. 

In applications of our method to measured cosmic rays data we therefore expect 
striking ehects if the galactic magnetic held parameterization and the directions of 
the sources are correct. 


6 Conclusions 

In this contribution we presented a method for evaluating dehections of cosmic rays 
in the galactic magnetic held. By making use of a held parameterization derived 
from Faraday rotation and synchrotron emission measurements we calculate the 
expected arrival directions of protons and use this for further analysis. 

To explore this method we introduced two observables designed for investigating 
cosmic ray dehections, namely an angular asymmetry and a measure of clustering 
strength. Applying the analysis to a simulated astrophysical scenario we demon¬ 
strated that the magnetic held corrections improve directional relations of cosmic 
rays and their sources if the correct held and source directions are used in the anal¬ 
ysis. This even works if the source directions are known with limited precision. 
If incorrect assumptions on the held or an arbitrary set of sources are used, the 
improvements in the observables remain marginal. In view of the galactic mag¬ 
netic held parameterizations based on measurements we thus expect that cosmic 
ray dehections in the galactic held can be verihed experimentally. 

In general, and independently of the AGN application, our method of magnetic 
held corrections enables studies of directional correlations between cosmic rays and 
messenger particles suitable to indicate directions of cosmic ray sources. 
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Figure 1: a) Simulated astrophysical scenario. The star symbols denote AGNs 
with the color corresponding to their distance. The circular symbols represent the 
cosmic rays (7% proton signal, 93% background contributions), b) Simulations of 
the expected arrival direction of protons after traversing the galactic magnetic held. 
The star symbol denotes the initial direction outside the galaxy, the cross symbols 
show the expected arrival directions of protons with different energies, and the color 
code gives relative probability distributions. 
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Figure 2: a) Angular distances of a cosmic ray to the nearest AGN direction and 
to the expected arrival direction assuming protons, b) Cumulated number of cosmic 
rays associated with AGN arrival directions with (symbols) and without magnetic 
held corrections (histogram) as a function of the maximum angular distance tto- c) 
Change in the angular distance before and after applying the magnetic held cor¬ 
rections. The vertical line separates reduced (> 0) and enlarged (< 0) angular 
distances. The vertical axis shows the smallest angular distance, whereas the hori¬ 
zontal line emphasizes the region below 5 deg. d) Frequencies of the cluster sizes m 
for the maximum angular distance 5 deg. 
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Figure 3: Improvement in the clustering strength versus the angular asymmetry 
using different analysis scenarios. The red cross denotes the data analysis using in 
a,c,d) the regular held, and in b) a reverse held. The box symbols denote the follow¬ 
ing variations: a,b) random held directions and isotropic cosmic rays, c) isotropic 
source directions, d) source directions varied by 15 deg. Refer to text for details. 
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